The coefficient of restitution does not exceed unity 

Hal Tasaki 2 

Abstract 

We study a classical mechanical problem in which a macroscopic ball is reflected 
by a non-deformable wall. The ball is modeled as a collection of classical parti- 
cles bound together by an arbitrary potential, and its internal degrees of freedom 
are initially set to be in thermal equilibrium. The wall is represented by an ar- 
bitrary potential which is translation invariant in two directions. We then prove 
that the final normal momentum can exceed the initial normal momentum at most 
by O(VmkT), where m is the total mass of the ball, k the Boltzmann constant, 
and T the temperature. This implies the well-known statement in the title in the 
macroscopic limit where O(VmkT) is negligible. Our result may be interpreted as 
a rigorous demonstration of the second law of thermodynamics in a system where a 
macroscopic dynamics and microscopic degrees of freedom are intrinsically coupled 
with each other. 



1 Introduction 

Impact between macroscopic bodies is a common phenomenon in which the complicated 
dynamics of microscopic degrees of freedom plays an essential role in determining the 
macroscopic motions of the bodies. A standard quantity characterizing impact is the 
coefficient of restitution eR, which is the ratio between the normal (relative) velocities 
before and after the impact. Although some elementary physics textbooks refer to the 
coefficient of restitution e R as a material constant, careful studies have revealed that eR 
depends sensitively on various conditions [IJ Ej . It has even been pointed out that the 
generally believed rule that eR cannot exceed unity may be violated in oblique impacts 
between a ball and a wall, provided that the wall is locally deformable jHlHI- 

Given such recent progress, it seems necessary to develop theories of impact from 
fundamental and statistical mechanical points of view. Here we prove a basic theorem 
which establishes the above mentioned rule eR < 1 (in the macroscopic limit) based on 
microscopic dynamics in a general and physically natural setting. More precisely we model 
a non-deformable wall by an arbitrary fixed potential which is translation invariant in two 
directions, and a ball as a collection of TV classical particles which are bound together by 
an arbitrary potential. We assume that the internal degrees of freedom of the ball are 
initially in their thermal equilibrium, and give the center of mass of the ball a macroscopic 
momentum. See Fig. Then we prove rigorous inequalities about the final momentum, 
which essentially implies eR < 1 for a macroscopic ball. Our bound contains a small 
correction (which is negligible for macroscopic bodies) which shows that even very small 
bodies can violate the eR < 1 rule only by a universal factor determined by the thermal 
momentum. 



^his version will appear in Journal of Statistical Physics. 
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Figure 1: A "ball" consisting of N classical particle is reflected by a non-deformable wall. 

The rule eR < 1 itself is by no means surprising. But as far as we know such a fun- 
damental and well known relation has not been established rigorously before. We believe 
that our explicit treatment of interplay between microscopic and macroscopic degrees of 
freedom in a dynamical system is unique and new. 

Of course the rule eR < 1 is closely related to the second law of thermodynamics. In 
the present problem (not in those treated in [HIE]), that eR < 1 essentially means that the 
internal thermal energy of the ball cannot be converted into macroscopic kinetic energy. 
This can be regarded as a manifestation of the second law in the form of Kelvin's law. We 
stress, however, that although there have been some rigorous derivations of the second law 
of thermodynamics based on microscopic mechanics El El HE CHI E], none of them 
apply to the present situation of a ball reflected by a wall. The reason for the inapplicability 
is essential rather than technical. In the existing proofs of the second law, macroscopic 
degrees of freedom are assumed to be completely controlled by an outside agent. The 
agent performs an operation (represented by a time-dependent Hamiltonian) according 
to a pre-fixed protocol, without allowing any feedback from the microscopic degrees of 
freedom 3 . In the present problem of impact, on the other hand, the macroscopic motion 
of the ball is determined only as a solution to the equation of motion involving all the 
degrees of freedom. It is therefore mandatory here to take into account feedback to the 
macroscopic motion from the microscopic motion. We accomplish this by making a full 
use of the symmetry of the model. 

Thus, from the view points of statistical mechanics and thermodynamics, our result 
eR < 1 can be regarded as a step in establishing the second law by properly taking into 
account feedback from microscopic degrees of freedom to macroscopic motion. Such a 
treatment is not only realistic but also sheds better light on the interplay between the 
macroscopic and the microscopic scales. We also note that in small physical systems (such 
as molecular motors) which operate under relatively large thermal fluctuation, feedback 
from internal degrees of freedom to overall motion plays an essential role. Recently Maes 
and the present author extended the present work, and proved similar (but weaker) results 
in a much more general setting fT^j . 

We note that the second law has been extended to nonequilibrium states in [T31 Ej 
by examining the Boltzmann entropy. Although these works do not treat concrete me- 

3 In thermodynamic settings, it is implicitly assumed that the feedback from microscopic degrees of 
freedom becomes negligible if the system and the means of operations become macroscopic. It is hard to 
judge how realistic such an idealization is. 
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chanical setup as our, they deal with very similar problem as the present work, namely, 
to demonstrate the second law by starting from a nonequilibrium state and following the 
deterministic dynamics. The constancy of the phase space volume plays an essential roles 
both in and the present work (as well as in the Boltzmann's original work). 

The reader might ask what our result imply about the thermalization process which 
should take place within the ball during and after the impact. Interestingly enough our 
analysis does not deal directly with possible thermalization, and still rigorously leads us 
to the desired rule < 1. This is indeed physically natural since a thermalization process 
is not universal while the rule eR < 1 is. Although it is likely that usually a sufficient 
thermalization takes place during the impact, this may not be the case when the impact 
occurs very quickly. Nevertheless the rule eR < 1 is always valid, and should be provable. 
It is of course challenging to study thermalization, but it requires a much more detailed 
and concrete analysis. 

2 Definitions and main results 

We consider a classical "ball" which consists of TV point-like particles whose canonical 
momenta and coordinates are denoted as and r^, respectively, for i — 1, 2, . . . , N. We 
use the Cartesian coordinates when expressing the momenta and coordinates in their 
components. The Hamiltonian of the ball is 

N \ \2 N 

h = E l P~ + £ u ^ - r *) + E v ^ (!) 

1=1 z / i * i>j 1=1 

where ^ is the mass of the i-th particle, Uij(-) is the potential of the force between the 
particles i and j, and V^(-) is the potential of the force that the wall exerts on the i- 
th particle. We assume that each V^(-) is invariant under any translation in the y or z 
directions. Other assumptions about the potentials of the wall will be stated later. As for 
the interaction, we assume that [/^(r^ — r^-) > Uo with a constant Uo, and Uijfa — rj) = oc 
when \ri — rj\ > Ro with a constant R > 0. The first assumption is introduced to 
guarantee that the canonical distribution (??) is well-defined. It may be replaced by a 
weaker condition if necessary. The second assumption guarantees that the ball cannot be 
larger than R . This constraint makes the proof less technical. Otherwise we have no 
requirements for Uij(-). 

We stress that the ball need not be a uniform one. One can design suitable Uij(-) (and 
constraints of the configuration if necessary) to treat, for example, an elastic sphere filled 
with a gas. Moreover the ball may have an arbitrary shape (determined by the potential 
Uij), so it need not be a "ball." 

As usual we introduce the momentum 4 

N 

P = (Px,P y ,Pz) := ^2pu (2) 
1=1 

4 We write A := B when we want to emphasize that B defines A. 
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and the coordinate 

1 N 

r = (x,y,z) := — J^r*, (3) 

of the center of mass, where m = X)* is the total mass of the ball. We also denote by 7 the 
collection of canonical coordinates for the remaining 3(N — 1) internal degrees of freedom. 
We denote these variables collectively as T = (p, r, 7) and denote by dT = dpdrd^f the 
Lebesgue measure on the whole 6iV-dimensional phase space. 

By using these new variables, we write the potential from the wall as 

N 

V(r, 1 ) = J2V i (r i ). (4) 

1=1 

Then the Hamiltonian (PJ is rewritten as 

H=^ + h^) + V(v n ), (5) 
2m 

where ^(7) is the Hamiltonian for the internal degrees of freedom. Let us list our assump- 
tions about the potential V(r, 7). 

i) V(x, y, z, 7) is independent of y and z. 

ii) V{x, y, z, 7) = if x < for any 7 such that ^(7) < 00. 

hi) For any 7 such that ^(7) < 00, one has j^V{x,y,z,^) > for any x > and 
-§^V(x, y ) z, 7) > fo for x > r . Here fo>0 and r > are constants. 

The translation invariance i) is a consequence of the already stated translation invari- 
ance of V^(r^). The property ii) means that the effect of the wall takes place only for 
sufficiently large x. As a consequence, the internal degrees of freedom of the ball and the 
motion of its center of mass are decoupled if x < 0. Then the center of mass moves with 
a constant velocity and the internal degrees of freedom 7 evolve separately according to 
the nontrivial internal dynamics determined by ^(7). The existence of such a "decoupled 
motion" plays an essential role in our proof. The property iii) means that the wall exerts 
repulsive force which becomes strong enough if the ball goes deeper into the wall. This 
guarantees that the ball is reflected by the wall. 

We denote by T(t) = (p(t), f (£), j(t)) a solution of the equation of motion determined 
by H. Of course this is a full fledged TV-body problem, and one can never write down the 
solution T(t) explicitly. Nevertheless we will be able to control some aspects of the motion 
rigorously. We denote the initial state as r (0) = (p(°), r^, 7^) := f(0). See Fig. [2(a). 

Let us specify the initial condition. Roughly speaking we take the initial state T^ ^ = 
(p(°) ? r(°), 7^) where the momentum p(°) of the center of mass (almost) has a macroscopic 
value po, and the internal degrees of freedom of the ball are in equilibrium. Therefore the 
ball initially has no macroscopic rotation. More precisely the initial state is sampled 
according to a measure p(T) dY ) where the weight p(Y) is written in a product form as 

p(r) = p m (p)p p (r) A ( 7 ). (6) 

Here the probability p p (r) for the position is an arbitrary normalized density (including 
delta functions) such that p p (x, y ) z) = if x > 0. Thus the ball does not feel the wall at 
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t — 0. As for the momentum we use for simplicity the normalized density given by 



n(Ti) - f3/{47r(Ap) 3 } if|p-p |<Ap; ( ) 

PmiPj 10 otherwise, [> 

where the initial macroscopic momentum is p = (po ,Pi,P2) with p > 2Ap > and 
arbitrary p±, p 2 - Finally we assume that the internal degrees of freedom of the ball are in 
equilibrium, and sample them according to the canonical distribution 

A(7) = W (8) 



with the partition function 



Z 



{(5) = Jd ie ~^\ (9) 



where (3 = (kT)~ l is the inverse temperature. Note that the canonical distribution is 
the physically correct choice (even for a very small ball) if the ball was in touch with a 
heat bath before the collision experiment 5 . When the ball can be regarded as a "healthy" 
macroscopic body, the choice of the equilibrium ensemble should not matter because of the 
"equivalence of ensemble." To demonstrate this, we treat in section El the situation where 
the internal degrees of freedom are initially distributed according to the microcanonical 
ensemble. 

Let if > be a time such that x(t) < for any t > if with probability one (for the 
initial state as described above). In other words the ball has certainly been reflected by 
the wall and is already in a "decoupled motion" at time if. Note that such if always 
exists because of the assumed property iii) about the force from the wall 6 . We shall fix 
if throughout the paper, and denote the final state as = (p^\r^\^) := T(if). See 
Fig. El (a). 

For any p > 0, let P{(p) be the probability that the final momentum = (p^\Py\pT ) ) 



satisfies \p%'\ > p. We define the thermal momentum as p t h = ym//3, which is extremely 
small for a macroscopic ball. Then the main result of the present paper is the following. 

Theorem: Set the uncertainty in the initial momentum as Ap = p^. Then one has 

p ' <p)si2l » !exp| -^ {p2 - to+p " )2}i - (io> 

and for the coefficient of restitution 

e R := < i + 2 ^ + 325 ( ^Vexp[-^]. (11) 

Po Po \PthJ Pth 



5 To be precise one must assume that the interaction between the ball and the bath is very weak, and 
that one can detach the ball from the bath without disturbing the state of the ball. 

6 We have assumed iii) only to ensure the existence of a finite if. Thus one can replace iii) with weaker 
conditions which still ensure the existence of if. 
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(p(l), r (l), 7 (l)) 



Figure 2: The solutions f (t) and f (t) are related by a shift of time and a translation. 
The variables in the figure play an important role in the proof. 

The bound f[TT|) shows that the coefficient of restitution e R does not exceed unity by 
much when p t ^ <C Po- Note that the condition p t ^ <C Po is always met in a macroscopic 
ball. The bound f[TUJ) contains stronger information, and shows that it is exponentially 
unlikely to have a final momentum p( f ) with \p$\ — Po ^ Pth- 

When the ball consists of a small number of particles and the thermal momentum p^ 
is not negligible (compared with p ), there can be some events 7 where the final momentum 
\p$\ becomes strictly larger than p . But our bound f[TU|l is valid even for such small balls, 
and establishes that the final momentum \p$ \ may exceed p G only by a factor comparable 
to pth- This means that when we say that a ball is macroscopic, we should speak about 
whether the initial momentum p is large compared with the thermal momentum p t \ V) 
rather than whether the ball consists of many particles. 

The requirement Ap = p^ in the theorem may sound natural from a practical point 
of view since it would be very hard to make Ap smaller than p^. Theoretically speaking, 
we expect that the inequalities (fTHj) and f[TT|) (with slight improvements) are valid in the 
limit Ap — > as well, but cannot prove them for purely technical reasons. The technical 
difficulty appears when we try to decouple macroscopic degrees of freedom (i.e., the p( f ) 
integral) from microscopic degrees of freedom (i.e., the 7^ integral) in (|23Jl . It is desirable 
to find an improved proof which works in the Ap —> limit. 

3 The case with a fixed protocol 

Before describing the proof, which is rather involved, let us describe the corresponding 
argument which works in the case where the macroscopic motion is treated as a fixed 
protocol as in the traditional proof of the second law [HI El El HE CHI E]- We hope that 
this short description will make clear both the essence of the proof and the difficulties 
encountered in the fully mechanical problem that we treat. 

Suppose that the internal degrees of freedom are governed by a time-dependent Hamil- 
tonian h(t\ 7) such that h(0] 7) = h(t{ ; 7) = ^(7). One may imagine that an outside agent 

7 One can easily construct such events with > Po, for example, in a model of two particles coupled 
by a harmonic potential. 
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performs a fixed operation on the ball. Again we denote by 7^ and 7^ the states at 
time t — and t — tf , respectively. From the energy conservation law, the energy that the 
external agent gains during the operation is 

AE :=h(^ 0) )-h(^). (12) 

Suppose that the initial state 7^ is distributed according to Pi(7^) of (JBJ). For any 
D > 0, let P{D) be the probability that AE > D. Then from the definition, one has 



e 



-/?M7 (0) ) 



P(D) = / rf 7 (0) - — - . (13) 
v ; Jae>d 1 Z((3) V ; 

Now the Liouville theorem implies d^^ = Then, from the energy conservation law 

(fT2j) . the probability f[T3j) can be bounded as 

r p -f3h(^))-f3AE r -ph(<yW)-pD 

p ^ - Lo^ ~^w^ £ l H 1 = e ~" D - <14) 

The final factor decays rapidly for D ^ 1//3, and implies that the energy gain AE cannot 
be larger than O(kT). In the macroscopic limit (where O(kT) is negligible), this implies 
Kelvin's version of second law, i.e., AE < 0. 

The above argument, which makes use only of the energy conservation law and the 
Liouville theorem, contains the basic strategy of our proof of the main theorem. But it is 
obvious that the same argument cannot be applied to our fully mechanical problem of the 
ball reflected by the wall. The main difficulty comes from the fact that the macroscopic 
motion itself becomes a dynamical degree of freedom, and the Liouville theorem does not 
imply d^ ^ = dj^\ but implies dT^ = dT^\ To overcome the difficulty, we shall make 
full use of the symmetry of the model and the property of the space of the momentum of 
the center of mass. 



4 Proof 

Note that we can set p± = P2 = by making a suitable Galilei transformation. Then the 
initial momentum is p = (po, 0, 0) plus a small fluctuation. 

Since the initial measure p(T) is normalized, the probability in question is written as 

Pt(p)= L, (0UI dT^pCrW), (15) 

J\p ( P(r(o))\> p 

where r( f )(rW) = (pW(r( c )),r( f '(r( '),7( f )(r( ))) is the final state represented as a function 
of the initial state . Since and are in one-to-one correspondence and the 
Liouville theorem implies dT^ = dT^\ one can rewrite ffTHj) as 

Pt(p)=[ mi ^ f ) p (r(°)(r( f ))), (16) 

j\p { P\>p 

where r^r*^) is the initial state represented as a function of the final state r^ f ^. Since 
r^ ) = (p(°), 7^) depends on = (p^,r^,7^) in a nontrivial manner, it is not 
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easy to evaluate the integral ffTHj) as it is. To "disentangle" the dependence as much as 
possible, we make a change of variables. 

Recall that T(t) = (p(t), ? (i), 7(*)) is a solution of the equation of motion with T(0) = 
r (0) and f (if) = T (f) . Since x (f) = x(t { ) < 0, there is a r > such that x(t f -r) = 0, where 
we wrote f (£) = (x(i), y(t), z(t)). Indeed there are two values of r satisfying x(tf — r) = 0, 
and we take the smaller one. In other words, the center of mass passes through the yz- 
plane (defined by x = 0) at t = U — r after the reflection from the wall. Since the motion 
"decouples" for t > U — r, the center of mass has a constant velocity /m, and we have 
r = mx^/pf. Let us write r (1) = (pW, r^, 7W) := f(t f - r). We then find pW = p( f ), 
r (i) = r (f) _ ( p (0/ m ) Tj and (by definition) x^ = 0. See Fig. [2 (a). 

Let us consider another solution T{t) = (p(£), f (£), j(t)) (which satisfies the same 
equation of motion determined by H) with the "final condition" T(tf) = (p^, o, 7^) where 
o = (0,0,0). We then trace back the time evolution and write = (p®,r®,7^) := 
f(0). See Fig. [21(b). We can naturally regard p®, y^ 2 \ and 7® as functions only of p^ 
and 7^. 

Note that the two solutions Y(t) and T{t) share the momentum p^ = p^ and the 
internal coordinate 7^ at the moment when the center of mass passes through the yz- 
plane after the impact. Since the whole system is translation invariant in the y and the 
z directions, this means that the two solutions are related by a simple shift of time and a 
spatial translation as p(t) = p(£ + r), f (£) = r(t + r) + r^, and j(t) = j(t + r). Noting 
that the motion "decouples" before the center of mass passes the yz-plane, these relations 
imply 

p(o)=p( ) = p(r)=p( 2 ), (17) 

and 

7 (0) =7(0) = 7(r) = F T ( 7 ( 2 )), (18) 

where F T (-) represents the time evolution generated by h(j) alone. As for the positions 
r(°) and r^ 2 \ we use r(r) = r(0) + (p/ra)r to get 

r (°) = r(0)=r(r)+r« 

D (0) 

= r(0) + — r + r« 

m 

(2)^P (0) + (f) 



m m 



,(2) 



+ r (f) + ^(p(o)_ p (f)). (19) 
m 



We then substitute r/m = /p$ to get 

r<°> = r( 2 ) + (0,</V f) ) + ^(pf ,pf -p£>,p?> -pf). (20) 

Suppose that p( f ) and 7^ are fixed, and consequently, p^ 2 ^ and r( 2 ) are also fixed. Then 
(|T5j> implies that the coordinate variables r^ ^ and r^ f ^ are related with each other by a 
simple linear transformation. This is the most important observation in the proof. 
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Since 7^ = ^(7^), the Liouville theorem for the internal degrees of freedom implies 
= d^ l \ We also have #(7^) = Pi (7®) since the canonical distribution is time- 
independent. By combining all these results we finally get 



Ptip) = f m dp^dr^d^p m (p {2) )p P (r)Pih {2) ) 

J \Px J \>P 



= [ (f) d P ^W 1] 



„(2) 
Px 



p m (p (2) )A(7 W ) 



(2)> 



(21) 



To get the second line, we substituted f|2Uj) and p(°) = p( 2 \ and integrated over r( f ) (with 
p( f ) and 7^ fixed). The factor \p$/p^\ is the Jacobian of the linear transformation (|20|) . 
Surprisingly, we were able to carry out the r^ f ^ integral exactly since p( 2 ), r^ 2 ^ and 7^ do 
not depend on r^ f ^ and p p (r) is normalized. The most delicate estimate is over. 
We now use the energy conservation law 



i P ( 2 )| 2 

2m 

and the explicit form of p\ to get 

Ptip) = L dp^d^ 1 

JpZ>>p 



+ /i( 7 (2) ) 



,( f )|2 



2m 



(22) 



< 



W\>p 
3 



V. 



(f) 



p. 



(2) 



P m ( P (2) ) 



-^(7 (1) )-^r{|P (f) l 2 -|P (2) l 2 } 



47r(Ap) 3 J\p(f)\> P 



dp (f) d 7 (1 



)_Wx 



e -/3^(7 (1) ) 



Po-Ap Z(0) 



-e 2 ™ 



lr{|p (f) l 2 -(po+A P ) 2 } 



(23) 



where we have enlarged the range of integral, and replaced p( 2 ) by its "worst possible 
values" and p m (p^) by its maximum value 3/{47r(Ap) 3 }. Now the p( f ) integral and the 
integral are completely separated and we get 



3 r |p( f) | 



Ap 



e -£{\p m \ 2 -( P0 +&p) 2 } 



Po{Ap) 6 



{p 2 -(po+Ap) 2 } 



(24) 



where we recalled that p t h = y m /ft, Po — ^P > Po/2, and bounded the integral using the 
standard (and elementary) estimate 



V IrVI „-|pT/(2a) - 



|p'|>P 



dp |p I e 



^'47r(p') 3 e- (p,)2/(2Q) 
ds4n(s + p) 3 e-^ 2 ^ 

< 47re- p2/(2a) I™ ds(s + pfe 

Jo 

< 8%y/a (p + a/«) 3 e 



3 _-« 2 /(2a) 



3 ^-p 2 /(2a) 



(25) 



where the final (crude) bound is checked by evaluating the integral explicitly. By setting 
Ap = pth in {2IJ) , we get the desired ([TO)) . 
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To deal with the coefficient of restitution e^, we note that for any S > 

< (|p (f) l) <po + s+ f rfr( f )|p( f )| P (r(r( f ))), (26) 



and bound the integral as before to get 



|p (f) |) < Po + 6 + 9 ^ ^(Po + ^ + Ah) 4 e -^ {(po+ , )2 _ (po+Art2}> (27) 



By setting 5 = 2Ap = 2p t h, we get (|TT|> . 



5 Microcanonical distribution 

We have proved our theorem in the setting where the internal degrees of freedom are ini- 
tially distributed according to the canonical distribution. To demonstrate that the choice 
of the equilibrium ensemble does not matter if the ball can be regarded as a "healthy" 
macroscopic system, we here show how to treat the case with the microcanonical distribu- 
tion. The following (not too short) analysis indeed corresponds to a single line calculation 
in (j2Sj) • From this one clearly sees that the canonical distribution is perfectly suited for a 
demonstration of the second law. 

Let the density of states of the internal degrees of freedom be 

<S>(E) := J d-fS(h(j)-E). (28) 

We assume that $(i?) satisfies the bounds 

d exp[Na(E/N)} < $(E) < c 2 exp[N a(E/N)], (29) 

with a function a(e) which is smooth, nonincreasing, and concave in e, and constants < 
C\ < c 2 . Here TV is the number of particles in the ball, and is assumed to be a large quantity. 
Of course a(e) is nothing but the entropy. We note that the above assumption may not hold 
in general. When this assumption is valid, the system exhibits "healthy" thermodynamic 
behaviors, and the canonical distribution and the microcanonical distribution become 
equivalent in the sense that they give the same thermodynamics in the limit TV — > oo. We 
shall now extend this equivalence to the case of the collision problem. 

We then choose two constants e > 0, AE > 0, and set the initial distribution of the 
internal degrees of freedom as 

Pi(l) = 7^ x[Ne < /i(7) < Ne + AE], (30) 

where the characteristic function is defined by x[true] = 1 and x [false] = 0. The normal- 
ization constant is 

rNeo+AE 

Q = / dE$(E) > a exp[7Va(e )] AE. (31) 

JNe 
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This defines the microcanonical distribution with a finite width AE. 

Now we substitute the present #(7) to the right-hand side of (|2T|) to get 



Pf(p) = 7^/ (f)| V f) rf7 (1) 
"0 -WI^p 

= ^/ (f) rf P (f) rf7 (1) 
"0 Wl^p 



pi 2) 



„(2) 



Pm(p (2) ) x[Ne < /i( 7 (2) ) < Ne + AE] 
p m (p (2) ) x 



X 



x[Ne < /i( 7 (1) ) + 



| p (f)|2 | p (2)|5 



2m 



2m 



<Ne + AE], (32) 



where we used the energy conservation (|22j> . We then note that for any p( 2 ) = p(°) such 
that p m (p^) 7^ 0, one has 



fa)) 2 
2m 



2m 



2m 



where 



A ^ = po Ap (Ap) 2 
m 2m 



(33) 



(34) 



Then © and (JHSD imply that 



^f(p) < 7^ / (f) rfp (f W X) 

"0 -WI^p 
x 



,(2) 



P m ( P (2) ) 



X 



x[iVe - A£ < %«) + - ^£ < 7Ve + AE + A£] 

zm zm 



p m ( P (2) ) 



•/ipi f) i>p 

x ^ / d7 (1) x[A^o - AE < %«) + - < iVeo + AE + AEp5) 

Note that the p^-* integral and the 7^) integral have been split. We can evaluate the 7^ 
integral separately as 



1 
1 



1) 



X [Ne - AE < h{^) + 



|p(f)|2 | p (o)p 



2m 



2m 



< Ne + AE + AE] 



< 



O 
1 



— / dE$(E) X [Ne — AE < E + 



| p (f)|2 | p (0)|2 



2m 



2m 



< Ne + AE + AE] 



1 



(AE + 2AE)c 2 exp[iV ( x(eo + -{ 
i Z iV zm 



p (f)|2 | p (0)| 



+ 



2m 



+ AE + AE})} 



< U + 



2AE\ c 2 exp[N a(e + N^{-^ + + AE + AE})] 



AE / 
2AE\ c 2 



AE J Cl 



— exp[/?{- 



ci exp[iV cr(e )] 

| p (f)|2 | p (0)|2 



2m 



+ 



2m 



+ AE + AE}], 



(36) 
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where we used (J3T]> to get the second inequality. To get the third inequality, we noted 
that (|2Sj) and the concavity of a(e) imply 

e N*(eo+(D/N)) < e Na(e )+f3D ^ /gyj 

for any L>, where /? = a'(eo). 

We now choose the uncertain AE (built into the microcanonical ensemble) comparable 
to AE ~ p Ap/m. Then by substituting (|5SJ| into (|S5j>. we get essentially the same 
estimate as in the first line of (|23j> with an extra constant. Note that the estimate gets 
worse as we take AE" smaller. 

It is a pleasure to thank Hisao Hayakawa, Tsuyoshi Mizuguchi, Hidetoshi Nishimori, 
Shin-ichi Sasa, and Ken Sekimoto for valuable discussions. I also thank referees for their 
constructive comments on the manuscript. 
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